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Abstract 


Observation of a large sample of HII galaxies shows that the emission line 
ratios of the youngest objects change systematically with gaseous oxygen abun- 
dance, which we interpret as resulting from changes in the initial mass function 
(IMF) of the ionising cluster. Comparison with cluster/nebula models shows that 
both the slope and the upper mass limit of the cluster IMF vary with abundance. 

In HII galaxies with oxygen abundance about 1/10 that of Orion, the IMF for mas- 
sive stars must have a slope which is about a factor of 2 smaller than in the So- 
lar Neighbourhood. 

Introduction 


"HII galaxies", dwarf galaxies with the spectra of giant HII regions, are 
ionised by massive 0 stars formed during a recent intense burst of star forma- 
tion. Although the very strong, narrow emission lines which characterise these 
objects have been used in several studies to determine physical conditions in 
the ionised region (eg Searle and Sargent 1972, French 1980, Kinman and David- 
son 1981, Kunth and Sargent 1983, Campbell, Terlevich and Melnick 1986a; CTMa) , 
little is known about their ionising clusters . 

The nebular excitation of low-abundance HII galaxies is higher than can be 
explained by a solar neighbourhood IMF (Lequeux et al 1981 , Bergvall 1985) . Gi- 
ant HII regions in spiral galaxies show trends in emission line ratios which in- 
dicate that the temperature of the ionising cluster, T tc>n , increases smoothly as 
the abundance decreases (eg Shields and Tinsley 1976, Stasinska 1980) . The ob- 
served change in T ton can arise if the IMF of the ionising cluster in giant HII 
regions and HII galaxies is abundance-sensitive: Shields and Tinsley (1976) 
proposed a systematic change in the upper mass limit of the IMF, while Terlevich 
(1982) suggested that the slope of the IMF decreases with decreasing abundance. 

Studies of the ionising cluster in HII galaxies have so far been based almost 
entirely on UV/optical/infrared colours, mass-to-light ratios, and recombina- 
tion line equivalent widths; the easily-observed emission line ratios, which 
carry important information about the mass distribution and age of the ionis- 
ing cluster, have been almost totally neglected. We report here the results of a 
comparison between the line ratios of a sample of ~ 50 HII galaxies and those of 
model nebulae ionised by clusters with a variable IMF. 

Modelling strategy 

The ionising clusters of HII galaxies evolve very rapidly, on timescales 
of order the lifetime of a massive star, a few Myr. As the cluster evolves, T ton 
falls and the "hardness" of the ionising spectrum decreases, leading to a rapid 
change in the relative strength of lines from different ionisation stages of a 
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given element. A least-squares fit to the distribution of HII galaxies in a T ton 
-sensitive emission line ratio will therefore always underestimate the T ton of 
the cluster IMF (as distinct from the cluster’s current T ton ) and consequently 
cause us to either overestimate the IMF slope or underestimate its upper mass 
limit. If we wish to study the abundance behaviour of the IMF, it is therefore 
crucial to remove the effect of cluster evolution. 

In principle, there exists a set of HII galaxies, covering a range in abun- 
dance, in which the ionising cluster is at close to zero age: a "zero age se- 
quence" (ZAS) which will appear as an upper envelope in T ton -sensitive emission 
line diagrams. Identification and modelling of this ZAS using zero-age main se- 
quence cluster models will allow us to determine the abundance behaviour of the 
HII galaxy IMF for stars more massive than ~ 20 M© . 

Model HII galaxies 

We have used the cluster code of Melnick, Terlevich and Eggleton (1985) to 
produce zero-age main sequence model clusters of mass 10 6 M© . The cluster mass 
distribution is specified by a power-law IMF with slope a between 0.1 M© and a 
maximum of 200 M© . We next pass the ionising part of the cluster’s spectrum to 
a model nebula code (Ferland and Truran 1981) and construct low-density, clumpy, 
radiation-bounded HII regions with oxygen abundances spanning the range covered 
by the data. The models are discussed in greater detail in Campbell, Terlevich 
and Melnick (1986b; CTMb) . 


Comparison of models and data: the ( fOIII) / [Oil] . 0/H) diagram 


The T ton -sensitive line ratio [0III]/[0II] is insensitive to 0/H and thus 
allows us to examine directly the behaviour of T ton with abundance. The ZAS, from 
which objects evolve downwards, is easily visible in figure 1. The curves in 
figure la are the loci of model HII galaxies ionised by clusters with a constant 
IMF specified by the M u and a shown. All these models are poor fits to the data. 
Those with a = 3.0 produce values of [0III]/[0II] which fall nearly an order 
of magnitude below the ZAS for 0/H a 4 x 10 -5 , while those with a - 1 .0 grossly 
overestimate the position of the ZAS for 6 x 10 -5 . It is immediately evident that 
for any M u , a "flat" IMF (a ~ 1.0) is required at low abundance. The steep turn- 
down of the [0III]/[0II] envelope with increasing abundance requires a system- 
atic decrease in M u or increase in a . 

The strongest dependence of a and M u on abundance which can be accommodated 
without conflicting with their observed values at the Solar Neighbourhood oxygen 
abundance are . „ 

M u = 5.9 (0/H) °' 3 (1) 

a = 1.2 log (0/H) + 6.9 (2) 


The loci of models in which M u and a vary according to equations (1) and (2) are 
shown in figures lb and lc respectively. The fit of the variable-M u models is 
little better than that of constant IMF models , owing to the small change in the 
cluster mass distribution caused by changes in M u . Variable a models do much 
better, but still overestimate the position of the envelope for 0/H ^1.5 x 10 -4 . 

It seems therefore that a systematic change with abundance of both a and M u 
is required to fit the [OIII] / [Oil] envelope. We have constructed two such mod- 
els, one in which a increases with abundance in the variable-M u model (equation 
(1)) , and another in which M u decreases with abundance in the variable-a model 
(equation (2)). Figure Id demonstrates that these closely similar models both 
produce a very good fit to the ZAS. 
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Figure 1 

The distribution of the sample in the ( [OIII] / [Oil] , 0/H) plane. Circles: data 
from CTMa; squares: data from Kunth and Sargent (1983) . IZwl8 is from Kinman and 
Davidson (1981) . The short dash-dot line represents the adopted zero-age enve- 
lope for 0/H > 1.5 x 10~ 4 , as explained in CTMb. The model loci are labelled as 
follows: C398 : constant IMF, a = 3.0, M u = 98 M© ; C3200 : constant IMF, a - 3.0, 
M u ■ 200 M© ; C198 : constant IMF, a = 1.0, M tt = 98 M© ; C1200 : constant IMF, a - 
1.0, M tt = 200 M© ; M3: IMF with a = 3.0 and variable M u (eq.l) ; Ml: IMF with a ■ 
1.0 and variable M u ; a 98 : IMF with M u « 98 M© and variable a (eq.2) ; a 200 : IMF 
with M u * 200 M© and variable a ; M a : the variable M u model (eq.l) with 1.0 < a 
< 2.7; a M : the variable a model (eq.2) with 200 > M u > 73 M© . 

Discussion 


An excellent fit to our adopted zero-age sequence in the ([OIII] /[Oil] , 

0/H) diagram has been obtained by allowing the mass distribution of the ionis- 
ing clusters of HII galaxies to vary systematically with abundance. Our models 
indicate a very pronounced change in the IMF in the oxygen abundance range 2.5 x 
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10 -5 to 2.5 x 10 -4 ; a is required to vary from c=: 1.5 to 2.7 and M u from 200 to 
a 75 M© . The dependence of the IMF slope on abundance is closely similar to that 
derived (using different methods) for giant HII regions and HII galaxies by Ter- 
levich (1982) . 

Conclusions 


Assuming trends in emission line ratios are due to changes in the stellar 
mass distribution of the ionising clusters of HII galaxies, the main results of 
the emission line modelling can be summarised as follows: 

1) No constant IMF can reproduce the values of the T ton -sensitive line ratio 
[0III]/[0II] observed in the youngest ("zero age”) objects over the oxygen abun- 
dance range 2.5 x 10 -5 < 0/H < 2.5 x 10 -4 . 

2) A systematic change with abundance of M u alone also fails to reproduce the 
zero age range in [0III]/[0II] , owing to the small effect which even large 
changes in M u have on the mass distribution of the ionising cluster. 

3) Models in which only a. varies cover much more of the zero-age range in 
[OIII] / COII] and if M u c* 200 M© , are correct for oxygen abundances Six 10 -4 . 
Above this they slightly overestimate the position of the [OIII] /[Oil] enve- 
lope . 

4) Models in which both a and M u vary with abundance reproduce the data envelope 
very well over the full range in oxygen abundance of 2.5 x 10 -6 to 2.5 x 10 -4 . The 
detailed abundance behaviour of oc and M u is given in CTMb. 

5) For the lowest -abundance HII galaxies, eg IZwl8 (0/H ~ 1.5x 10 -5 ) , the pre- 
dicted IMF slope is o* 1.2, irrespective of M u . IZwl8 is an evolved object in 
which T ton is currently considerably lower than its zero-age value. 
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